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We analyze the electromagnetic form factors of the nucleon to fourth order in rela- 
tivistic baryon chiral perturbation theory. We employ the recently proposed infrared 
Mh. regularization scheme and show that the convergence of the chiral expansion is im- 

(— I I proved as compared to the heavy fermion approach. We also discuss the inclusion of 

vector mesons and obtain an accurate description of all four nucleon form factors for 
, . ■ momentum transfer squared up to Q^ ~ 0.4 GeV^. 
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1 Introduction 

The electromagnetic structure of the nucleon as revealed in elastic electron-nucleon scat- 
tering is parameterized in terms of four form factors]^ The understanding of these form 



factors is of utmost importance in any theory or model of the strong interactions. Abundant 
data on these form factors over a large range of momentum transfer already exist, and this 
data base will considerably improve in the few GeV region as soon as further experiments 
at CEBAF will be completed and analyzed. In addition, experiments involving polarized 
beams and/or targets are also performed at lower energies to give better data in particular 
for the electric form factor of the neutron, but also for the magnetic proton and neutron 
ones. Such kinds of experiments have been performed or are under way at NIKHEF, MAMI, 
ELSA, MIT-Bates and other places. Clearly, theory has to provide a tool to interpret these 
data in a model-independent fashion. For small momentum transfer, this can be done in the 
framework of baryon chiral perturbation theory (ChPT), which is the effective field theory of 
the Standard Model at low energies. This will be the main topic of the present investigation. 

To put our work presented here into a better perspective, let us recall what is already known 
from chiral perturbation theory studies of the electromagnetic form factors. Already a long 
time ago it was established that the isovector (Dirac and Pauli) charge radii diverge in the 
chiral limit of vanishing pion mass [0]. The first systematic investigation in the framework 
of relativistic baryon ChPT was given in p|, in particular, it was shown that the analytic 
structure of the one-loop representation of the isovector spectral functions is in agreement 
with the one deduced from unitarity (in the low energy region, that is on the left wing of the 
rho resonance). That approach, however, suffers from the fact that due to the use of standard 
dimensional regularization, the one-to-one correspondence between the expansion in loops 
and the one in small momenta was upset. If one considers the nucleon as a very heavy 
static source, a consistent power counting is possible, the so-called heavy baryon chiral 
perturbation theory (HBChPT). Within that approach, the electromagnetic form factors 
were studied in 0, ^, the latter also containing the extension to an effective field theory 
including the delta resonance. It was found that the chiral description already fails for values 
of the four-momentum transfer squared of about Q^ ~ 0.2 GeV^. In 0, the isovector and 
isoscalar spectral functions were investigated in the heavy nucleon approach. It was pointed 
out that due to the heavy mass expansion, the analytic structure of the isovector spectral 
function is distorted, making the chiral expansion fail to converge in certain regions of 
small momentum transfer. This can be overcome in the recently proposed Lorentz-invariant 
formulation of making use of the so-called "infrared regularization". Being relativistic, 
this approach leads by construction to the correct behavior of the spectral functions in the 
low energy domain. Furthermore, it is expected to improve the convergence of the chiral 
expansion. This will be one of the main issues to be addressed here. We perform a complete 
one-loop analysis of the form factors, i.e. taking into account all terms up to fourth order. 
We will demonstrate that one can achieve a good description of the neutron charge form 



#''As will be discussed in more detail later, one can either work with the Dirac and Pauli form factors Fi 
and F2 or the Sachs form factors Ge and Gm of the proton and the neutron. 



factor for momentum transfer squared up to about Q^ = 0.4 GeV^. The other three form 
factors cannot be precisely reproduced by the pion cloud plus local contact terms to this 
order. The source of this deficiency is readily located - it stems from the contribution of 
vector mesons. We show how to incorporate these in a chirally symmetric manner without 
introducing additional free parameters (the new parameters related to the vector mesons 
are taken from dispersion-theoretical analyses of the form factors, see [0, ^ 0). Within 
that framework, we obtain a very precise description of the large "dipole-like" form factors 
without destroying the good result for the neutron form factor obtained in the pure chiral 
expansion. It is also important to stress that the results obtained in the heavy baryon 
approach can be straightforwardly deduced from the framework employed here, shedding 
some more light on previous HBChPT results. 

The manuscript is organized as follows. In section 0, some basic definitions concerning the 
electromagnetic form factors are collected. The one-loop representation of the nucleon form 
factors is given in section ^ The effective Lagrangian on which our investigation is based is 
briefly reviewed in subsection p.l| . In particular, infrared regularization and the treatment 



of loop integrals are discussed in some detail in subsection p^ The pertinent results based 
on this complete one-loop representation are presented and discussed in subsections ^]3] and 
|3.4| . Particular emphasis is put on a direct comparison with the results obtained in heavy 
baryon chiral perturbation theory, which is a limiting case of the procedure employed in this 
manuscript. The inclusion of vector mesons in harmony with chiral symmetry is presented in 
section ^. We give some technicalities in subsection [4.1| and display the results for the form 
factors in E^, followed by some remarks on resonance saturation in [4.3|. Concluding remarks 



are given in section ^. Some further technical aspects are relegated to the appendices. 

2 Nucleon form factors 

The structure of the nucleon (denoted by W) as probed by virtual photons is parameterized 
in terms of four form factors, 

{Nip') I J, I Nip)) = euip') [i.F^it) + '.^F^it)Yip) , N = p,n, (2.1) 

with t = q^q^ = (p' — p)"^ the invariant momentum transfer squared, J7^ the isovector vector 
quark current, J'^ = qQ'J^q (Q is the quark charge matrix), and m^v the mean nucleon mass. 
In electron scattering, t is negative and it is often convenient to define the positive quantity 
Q^ = —t > 0. Fi and F2 are called the Dirac and the Pauli form factor, respectively, with 
the normalizations -Ff (0) = 1, -F"(0) = 0, -FflO) = ^p ^^^ -^2^(0) = '^n- Here, k, denotes the 
anomalous magnetic moment. One also uses the electric and magnetic Sachs form factors. 
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GEit) = Frit) + ^3^F2(t) , GMit) = Frit) + F^it) . (2.2) 
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In the Breit-frame, Ge and Gm are nothing but the Fourier-transforms of the charge and the 
magnetization distribution, respectively. In a relativistic framework, as it is used throughout 



this text, the Dirac and Pauh form factors arise if one constructs the most general nucleonic 
matrix element of the electromagnetic current consistent with Lorentz invariance, parity and 
charge conjugation. In the non-relativistic limit on the other hand, in which the nucleon 
can be considered as a very heavy static source, one naturally deals with the Sachs form 
factors. Therefore, as first stressed in f^, these arise in the heavy baryon approach to chiral 
perturbation theory. Since we will frequently compare our results to those obtained in that 
framework, we will consider both Pauli and Dirac and the Sachs form factors. For the 
theoretical analysis, it is advantageous to work in the isospin basis, 

Frit) = Fm±FPit), (z = l,2), (2.3) 

since the photon has an isoscalar (J = s) and an isovector (/ = v) component (and similarly 
for the Sachs form factors). It is important to note that the lowest hadronic states to which 
the isoscalar and isovector photons couple are the two and three pion systems, respectively. 
These are the corresponding thresholds for the absorptive parts of the isoscalar and isovector 
nucleon form factors. 

The slope of the form factors at t = is conventionally expressed in terms of a nucleon 
radius (r^)^/^, 

F(t)=F(0)(l + l(r2)t + ...) (2.4) 

which is rooted in the non-relativistic description of the scattering process in which a point- 
like charged particle interacts with a given charge distribution p(r). The mean square radius 
of this charge distribution is given by 

6 dF{Q^) 



{r ) = Att drr p{r) 



(2.5) 



/o ^' ' F(0) dQ^ 

Eq. ( |2.5|) can be used for all form factors except G% and F^ which vanish at t = 0. In these 
cases, one simply drops the normalization factor 1/F(0) and defines e.g. the neutron charge 
radius via 

dG%m 



{(rlY) =-Q 



dQ^ 
It is important to note that the slopes of G^ and F" are related via 

dGUQ') dF^Q') F-iO) 



(2.6) 



dQ 
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2=0 dQ^ 

where the second term in eq. (|2.7| ) is called the Foldy term. It gives the dominant contribution 
to the slope of G%. 

The large body of electron scattering data spanning momentum transfers from Q^ ^ to 
Q"^ ^ 35 GeV^ can be analyzed in a largely model-independent fashion using dispersion 
relations 0, |, 0. Here, we are interested in the region of small momentum transfer, Q"^ < 
0.5 GeV^. Nevertheless, since there is still some substantial scatter in the data in this range 
of momentum transfer, we will also use the results of the dispersive analysis for comparison 
to the ones obtained in the chiral expansion. A recent review on the theory of the form 
factors is given in [0, the status of the data as of 1999 is discussed in |jTll . 



3 One— loop representation 

3.1 Effective Lagrangian 

Our starting point is an effective field theory of asymptotically observable relativistic spin- 
1/2 fields, the nucleons, chirally coupled to the pseudo-Goldstone bosons of QCD, the pions, 
and external sources (like e.g. the photon field). The theory shares the symmetries of QCD 
(spontaneously and explicitly broken chiral symmetry, parity, charge conjugation, time re- 
versal invariance, and Lorentz covariance) and can be formulated in such a way as to obey 
systematic power counting, which will be discussed at length in the following section. Exter- 
nal momenta and quark (pion) mass insertions are treated as small quantities in comparison 
to the scale of chiral symmetry breaking, A^ ~ 1 GeV. The most direct link of the effective 
field theory to the underlying one, QCD, comes from the analysis of the chiral Ward iden- 
tities as stressed by Leutwyler |T^. Any contribution to S-matrix elements or transition 
currents has the form 

-M = g^/(f,^) , (3.1) 

where g is a generic symbol for any small quantity, / a function of order one, A a regularization 
scale, and g a collection of coupling constants. Chiral symmetry demands that the power i/ 
is bounded from below. The expansion in increasing powers of this parameter is called the 
chiral expansion. At a given order, one has to consider tree as well as loop diagrams, the 
latter ones restoring unitarity in a perturbative fashion. The machinery to perform this is 
based on an effective Lagrangian, which consists of a string of terms of increasing (chiral) 

dimension, 

r - r(2) I r(i) _l r^^) , r^^) , rW , /o on 

where the ellipsis denotes terms of higher order not needed here. We note that we perform a 
complete one-loop analysis, i.e. taking tree graphs with insertions from all terms indicated 
in eq. ( |3.2| ) and one-loop diagrams with insertions from £^^ and at most one insertion from 

We will now consider the terms relevant to our problem (for a more detailed description, 
we refer the reader to [0). The chiral effective pion Lagrangian, which to leading order 
contains two parameters, the pion decay constant (in the chiral limit) F and the pion mass 
(its leading term in the quark mass expansion) M, is given by 

4'i = ^(^X + X+), (3.3) 

where the triplet of pion fields is collected in the SU(2) valued matrix U{x) = u^(x), and 
the chiral vielbein is related to u via m^ = i{u\ V^m}. V^ is the covariant derivative on the 
pion fields including external vector (t>^) and axial (a^) sources, V^^ = d^U — i{Vf^ + af^)U + 
iU{v^ — Qfj,). The mass term is included in the field x+ via the definitions x = '^B{s + ip) 
and x± = u^X'>J'^ =t ''J'X^'^y with s and p being scalar and pseudoscalar sources, respectively, 
the former including the quark mass matrix, s = Ai + . . . . Furthermore, (. . .) denotes the 
trace in fiavor space. 



The pion-nucleon Lagrangian at leading order reads 

/:a = ^(z^-m + ^^75)^, (3.4) 

where the bi-spinor \l/ collects the proton and neutron fields and qa is the axial-vector 
coupling constant measured in neutron /3-decay, qa = 1.26. To be more precise, the nucleon 
mass and the axial coupling should be taken at their values in the two fiavor chiral limit 
{rriu = rrid = 0, rUs fixed at its physical value). To this order, the photon field only couples 
to the charge of the nucleon. It resides in the chiral covariant derivative, -D^\l/ = 9^\1/ + r^\l/, 
with the chiral connection given by F^ = \[u\ 9^-u] — |-u"''(t>^ + a^)M — ^u{v^ — a^)u^ . 

The minimal Lagrangians at second and third order have been given in |T^. We only show 
the terms needed for the calculation of the form factors, 



m/.,M^ 



4'il = ^ {ci(x+) - ^,{{u,uJ){D^, D-^} + h.c.) + ^a^^[t 

+ ^^"'F^. + ^^""{K^)] ^ , (3-5) 

8m ^ 6171 ^ J 

4'il = ^{^([/^^F+p. + h.c.) + ^([D^(F+)p. + h.c.)}vl/. (3.6) 

Here, F^ = u^F^yU + uF^^u\ and F^,^ = d^Ai, — dyA^ is the conventional photon field 
strength tensor. Furthermore, we have adopted the notation introduced in |T^ for traceless 
operators in SU(2), A = A — ^{A). The Cj and di are the so-called low energy constants 
(LECs) , which encode information about the more massive states not contained in the effec- 
tive field theory or other short distance effects. These parameters have to be pinned down 
from some data. In our case, the LECs ci, C2, C4, which only appear in loop diagrams, can 
be taken from analyses of vrA^-scattering, whereas ce, cj parameterize the leading magnetic 
photon coupling to the nucleon, and de, dj have to be fitted to the charge radii of proton 
and neutron. 



The minimal pion-nucleon Lagrangian to fourth order has been worked out in |]T5| . Of the 
118 terms given there, we only show the four which are of interest here. 



ay 



-'-f{K^) (X+) ^'"^ - ^^i(x+) ^''] ^ • (3.7) 

Two more terms (numbered 107 and 108 respectively in |]I5[) only contribute when taking 
into account effects due to m^, 7^ nid- We shall disregard these in what follows. We also 
note that the terms ~ 6105, Cioe only amount to a quark mass renormalization of the leading 
magnetic couplings ~ Cg, C7, 

C6 -^ ce = ce - 16 m M^ cioe , 

cy -^ C7 = C7-8mM^(2eio5 -eioe) . (3.8) 

When it comes to numerical evaluation, we will just use the renormalized constants cq and 
£7 and will not regard eios and eioe as additional parameters to be fitted. 
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3.2 Infrared regular izat ion 

Chiral perturbation theory in the meson sector as an effective theory for weakly interacting 
Goldstone bosons allows for a systematic expansion of physical observables simultaneously 
in powers of small momenta and quark masses. For example, in the isospin limit {rriu = m^) 
the quark mass expansion of the pion mass takes the form 

^' = ^'|l - ^^^^^^"4 + 0{M^) , (3.9) 

where the renormalized coupling £3 depends logarithmically on the quark mass, M'^di^/dM^ 
= —1. The infinite contribution of the pion tadpole ~ ^/{d — 4), see graph (a) in fig. |l|, has 
been absorbed in the infinite part of this LEG. In this scheme, there is a consistent power 
counting since the only mass scale (i.e. the pion mass) vanishes in the chiral limit rriu = 
rrid = 0. If one uses the same method in the nucleon case, one encounters a scale problem: 
for instance, if the nucleon field is treated relativistically based on standard dimensional 
regularization, the nucleon mass shift calculated from the self-energy diagram (see graph 
(b) in fig. |I|) can be expressed via [Q] 

mN-m= ^ m I Co + Ci/i^ - TCfi^ - /i^ log /i + ^ a^fi" | , (3.10) 

with /i = M/m, and Cq, Ci are renormalized LECs (the precise relation of which to previously 
defined LECs is of no relevance here). This expansion is very different from the one for the 
pion mass, in that the nucleon mass already receives a (infinite) renormalization in the chiral 
limit. This difference is due to the fact that the nucleon mass does not vanish in the chiral 
limit and thus introduces a new mass scale apart from the one set by the quark masses. 
Therefore, any power of the quark masses can be generated by chiral loops in the nucleon 
case, whereas in the meson case a loop order corresponds to a definite number of quark 
mass insertions. This is the reason why one has introduced the heavy mass expansion in the 
nucleon case. Since in that formalism the nucleon mass is transformed from the propagator 
into a string of vertices with increasing powers of l/m, a consistent power counting can be 
formulated. However, this method has the disadvantage that certain types of diagrams are 
at odds with strictures from analyticity. The best example is the so-called triangle graph, 
which enters e.g. the scalar form factor or the isovector electromagnetic form factors of the 
nucleon (see graph (6) in fig. ^. This diagram has its threshold at to = 4M^ but also a 
singularity on the second Riemann sheet, at tc = 4:M^ — M^/mj^ = 3.98M^, i.e. very close to 
threshold. To leading order in the heavy baryon approach, this singularity coalesces with the 
threshold and thus causes problems (a more detailed discussion can be found e.g. in p, p!6|). 
In a fully relativistic treatment, such constraints from analyticity are automatically fulfilled. 



It was recently argued in |]T^] that relativistic one-loop integrals can be separated into "soft" 
and "hard" parts. While for the former, a similar power counting as in HBChPT applies, the 
contributions from the latter can be absorbed in certain LECs. In this way, one can combine 
the advantages of both methods. A more formal and rigorous implementation of such a 



program is due to Becher and Leutwyler [^ . They call their method, which we will also use 
here, "infrared regularization" . Any one-loop integral H is split into an infrared singular 
and a regular part by a particular choice of Feynman parameterization. Consider first the 
regular part, called R. If one chirally expands these terms, one generates polynomials in 
momenta and quark masses. Consequently, to any order, R can be absorbed in the LECs 
of the effective Lagrangian. On the other hand, the infrared singular part / has the same 
analytical properties as the full integral H in the low-energy region and its chiral expansion 
leads to the non-trivial momentum and quark-mass dependences of ChPT, like e.g. the chiral 
logarithms or fractional powers of the quark masses. It is this infrared singular part I that is 
closely related to the heavy baryon expansion: there, the relativistic nucleon propagator is 
replaced by a heavy baryon propagator plus a series of 1/m-suppressed insertions. Summing 
up all heavy baryon diagrams with all internal-line insertions yields the infrared singular 
part of the corresponding relativistic diagram. 

To be specific, consider the self-energy diagram (b) of fig. |l]. In d dimensions, the corre- 
sponding scalar loop integral is 

^^^^^~i I {2txY [A'P -k^- te] [m2 - (p - kf - ie] ' ^^'^^'^ 

At threshold, p^ = sq = (M + m)^, this results in 

H{so)=c{d) ^ + ^ =/ + i?, (3.12) 

with c{d) some constant depending on the dimensionality of space-time. The infrared sin- 
gular piece I is characterized by fractional powers in the pion mass and generated by loop 
momenta of order M^^. For these soft contributions, the power counting is fine. On the other 
hand, the infrared regular part R is characterized by integer powers in the pion mass and 
generated by internal momenta of the order of the nucleon mass (the large mass scale). These 
are the terms which lead to the violation of the power counting in the standard dimensional 
regularization discussed above. For the self-energy integral, this splitting can be achieved 
in the following way: 



{2tiYAB Jo J {2nY[{l-z)A + zB] 



2 



^^/(tn(l-z)A + zi?P=' + ^' ^'•'') 

with A = M^ — k"^ — ie, B = m? — (p — kY — ie. Any general one-loop diagram with 
arbitrary many insertions from external sources can be brought into this form by combining 
the propagators to a single pion- and a single nucleon-propagator. It was also shown that this 
procedure leads to a unique, i.e. process-independent result, in accordance with the chiral 
Ward identities of QCD. This is essentially based on the fact that terms with fractional versus 
integer powers in the pion mass must be separately chirally symmetric. Consequently, the 
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transition from any one-loop graph H to its infrared singular piece / defines a sjTiinietry- 
preserving regularization. However, at present it is not known how to generalize this method 
to higher loop orders. Also, its phenomenological consequences have not been explored in 
great detail so far. It is, however, expected that this approach will be applicable in a larger 
energy range than the heavy baryon approach. 

Some remarks concerning renormalization within this scheme are in order. To leading order, 
the infrared singular parts coincide with the heavy baryon expansion, in particular the infi- 
nite parts of loop integrals are the same. Therefore, the /3-functions for low energy constants 
which absorb these infinities are identical. However, infrared singular parts of relativistic 
loop integrals also contain infinite parts which are suppressed by powers of /i, which hence 
cannot be absorbed as long as one only introduces counter terms to a finite order: exact 
renormalization only works up to the order at which one works, higher order divergences 
have to be removed by hand. Closely related to this problem is the one of the new mass 
scale A which one has to introduce in the process of regularization and renormalization. In 
dimensional regularization and related schemes, loop diagrams depend logarithmically on 
A. This log A dependence is compensated for by running coupling constants, the running 
behavior being determined by the corresponding /3-functions. In the same way as the con- 
tact terms cannot consistently absorb higher order divergences, their /3-functions cannot 
compensate for scale dependence which is suppressed by powers of /i. In order to avoid 
this unphysical scale dependence in physical results, the authors of [P| have argued that the 
nucleon mass mj^ serves as a "natural" scale in a relativistic baryon ChPT loop calculation 
and that therefore one should set A = ttin everywhere when using the infrared regularization 
scheme. This was already suggested in ^ for the framework of a relativistic theory with 
ordinary dimensional regularization. We will follow this idea, hence in our formulae, what 
would actually be log(M,r/A) will always be spelled out as log/i. 

In the following calculation, we frequently compare our results to the equivalent heavy baryon 
ones. We wish to emphasize that one can always regain the heavy baryon result from the 
infrared regularized relativistic one by performing a strict chiral expansion of all involved 
loop functions. As a simple example, we again show the expansion of the nucleon mass up 
to third order corresponding to eq. ( p. 101) , but this time using infrared regularization. The 



result is 

m.T^T — m. ^ — dr-, A/f 4- 

2F2 



mj^-m = -4ciM2 + JiAmM^I{m'^) , (3.14) 



where the loop integral I{ra ) is given by 



A-') = -/^'(;^ + T^l°g/^j + T^|f-\/4-/^^arccos(^-|jj. (3.15) 

(Definitions and results for this and all other loop functions needed in this work can be found 
in appendix ^.) Expanding to leading order, one finds 

/(m^) = --^ + 0(/i2) , (3.16) 



which leads to the well-known result 



?,g\M^ , ,,, 4^ 



m^-m = -Ac,M' - -^^ + 0{fi^) , (3.17) 

yielding a contribution non-analytic in the quark masses. As explained before, the loop 
function livr?) contains a non-leading divergence which cannot be absorbed to this order, 
but will be by an appropriate contact term at fourth order. The calculation of the nucleon 
mass to fourth order and the reduction to the heavy baryon limit is done in detail in 0. 

An essential ingredient to the treatment of loop integrals as described in the previous para- 
graphs is the fact that higher order effects are included as compared to the "strict" chiral 
expansion in the heavy baryon formalism. This was justified in |^ by improved convergence 
properties in the low energy region, it introduces however a certain amount of arbitrariness 
as to which of these higher order terms to keep and which to dismiss. It is therefore manda- 
tory to exactly describe our treatment of these terms. The philosophy in 0] was, above 
all, to preserve the correct relativistic analyticity properties. This was achieved by keeping 
the full denominators of loop integrals (and evaluating them by the infrared regularization 
prescription), while expanding the numerators to the desired chiral order only. In addition, 
e.g. crossing symmetry is to be conserved. We explore a different approach here: as we 
anticipate the neutron electric form factor to be highly sensitive to recoil effects, we keep 
all terms which occur according to the infrared regularization prescription and do not even 
expand the numerators of loop integrals. However, three effects involving low energy con- 
stants have to be discussed separately, all in the spirit of avoiding "artificial" higher order 
counterterms: 

• Loop diagrams with second order insertions proportional to c\ can easily be summarized 
by a shift of the "bare" nucleon mass to its renormalized value at second order, m — »• 
mAT = m — 4ciM^ + 0[(^\ In principle, this does not occur in those loop diagrams 
which have insertions from other second order low energy constants. As we do not 
want to artificially introduce this additional LEG by discriminating between m and 
nifq in our formulae, we allow for this renormalization everywhere. 

• As detailed in the previous section, the second order LECs ce, cy receive a renormal- 
ization ~ M^ at fourth order. In principle, to this order only the unrenormalized Cj 
appear in loop diagrams. We disregard this difference and use the same values both 
on tree level and for the loops. 

• From the definitions of the various electromagnetic form factors, it is obvious that these 
are not all calculated to the same accuracy in terms of chiral orders. A calculation 
employing chiral Lagrangians up to and including 0(5^) will be able to give the Dirac 
form factor F^ to 0{(f'^ and the Pauli form factor F^ to 0{(^\ When combining these to 
the Sachs form factors Ge and Gm-, the chiral orders are "mixed", see eq. (p^). When 
we, in the following, talk about "third" / "fourth" order calculations of the various 
form factors, we always refer to the projections of the third/fourth order amplitude 
onto these observables. We truncate exactly one kind of term in this process: the 
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0{q^) counterterms which enter the electric (charge) radii are really "G^; counterterms" 
in the sense that they appear in Fi and F2 with opposite signs in order to cancel 
exactly in Gm- In addition, there are, at C(g^), F2 counterterms fixing the magnetic 
radii. Therefore, Ge receives g'* counterterms inherited from (g^/4rri^)F2. As a genuine 
polynomial contribution of this kind only appears in an 0{q^) calculation, we drop 
these terms in Ge- 



3.3 Chiral expansion of the nucleon form factors 

The chiral expansion of a form factor F (being a genuine symbol for any of the four elec- 
tromagnetic nucleon form factors) consists of two contributions, tree and loop graphs. The 
tree graphs comprise that of lowest order with a fixed coupling (the nucleon charge) as well 
as counterterms from the second, third, and fourth order Lagrangians. As one-loop graphs, 
we have both those with just lowest order couplings and those with exactly one insertion 
from C\pf. The pertinent tree and loop graphs are depicted in fig. ^ (we have not shown 
the diagrams leading to wave function renormalization) . The first four of these comprise the 
tree graphs with insertions up to dimension four, and diagrams (5) to (9) ((10) to (12)) are 
the third (fourth) order loop graphs. Consequently, any form factor can be written as 



F{t) = F^'^'^t) + F^°°P{t) 



(3.18) 



In appendix P, we have listed the contributions to the Dirac and Pauli form factors from the 
various diagrams, including the nucleon Z-factor. This allows to reconstruct the pertinent 
expressions for Fi'2 in a straightforward manner. We refrain from giving the complete 



expressions here. The corresponding Sachs form factors G^V are obtained by use of eq. 
under the restrictions discussed in the preceding section. We only show the explicit formulae 
for the magnetic moments, the electric and the magnetic radii (in terms of renormalized, i.e. 
physical, masses and coupling constants). The magnetic moments read 



K 



+ 



(47rF, 



mr'<^-i) 



ce 



+ 



arccos 



M^ r 



AiV4 - /i^ 
2(7-3^^2 



2(^8-13/i^ + 3/i^ 
5 



/^' ( 8 - ^/^' I C6 



4 - -/i' ) C6 



log/i 



ji^mlf ( 3 



I r 3 't 

J^l -l/rriN C2 - (Sf/mN C2 - Stun C4 + 2cej log/x | , (3.19) 



(47rF, 
cl + 2c^7 

(AnF^y |Yi2 + C6 + 2c, 



yuarccosi^— ^ 
V4-/x^ 



M^ r 



2(3-^2^ + (4-^^2)(c6 + 2c7 
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+ 



3/i^m^ C2 / 1 , 

log/i 



log/i 



(47rF^)2 V4 



(3.20) 



with fi = MTr/rrijsf ^ 1/7. Note that the magnetic moments are finite in the chiral hmit. We 
remind the reader of the definitions of Cg and cy in eq. ( |3.8|) . We have distinguished between 
the renormahzed and the unrenormahzed LECs (at tree level and in loop contributions, 
respectively) only for reasons of formal exactitude, and will disregard this difference from 
now on. Note however that, in contrast to cq, cy, these renormahzed LECs inherit infinite 
parts from eios, eioe which are needed in order to remove the leading order divergences in 
the magnetic moments (as explained in section |3.2| ). Their poles in {d — 4) are determined 
by the corresponding /3-functions, 



Peios 
Peioe 



3^i 



16^. (^ + ^^ + 2^0' 
\7g\ - ArriNCi + (1 + 2g\)c^ 



irriN 



(3.21) 
(3.22) 



Cg, Cy in eqs. ( |3.19| ), ( |3.2(]| ) denote the finite LECs, i.e. with the poles subtracted. Next, the 
charge radii can be brought into the form 



^v \2 



-I2d' + 



9l 



2m% 



ivrF^ 



ifi 



4^^(1 



C6 



yUarccos(— ^^ 
v/4^7? 



70 - 74^^ + 15;u^ + 3/i^ 3 - fi')ce 



„s \2 



+ 



(47rF^) 
-24 dj + 



10 - 44/^2 + 15/i^ 



l + 21og/i 



3/i^(l - fi'^jCQ 



log/i 



(3.23) 



(47rF^) 



3 k' 
2m% 



/i 



^ 4 



1 + 3 C6 + 2c, 



yuarccosi^— ^ 

(4 - /i2)3/2 



54 - 34/^2 + 5/i^ 



+ 



4-5/2^ + 3 l-/i 



(c6 + 2c^ 



-3(l2-7/i2 + /i')(c6 + 2c, 

logyU 



(3.24) 



The isovector charge radius diverges logarithmically in the chiral limit, while the isoscalar one 
is finite, d^ denotes the renormahzed LEC t/g, i.e. without the infinite part that removes the 
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leading order divergence in the isovector electric form factor, the corresponding /3-function 
being given by [|18|, |14 

/?., = -^^ . (3.25) 

The LEG d-j entering the isoscalar radius is, in contrast, finite. Finally, the magnetic radii 
are given by 



^^''{rMf 



74 



+ 



24mAre 



+ 



arccosi^— ^j 

Ai(4 - /i2)3/2 



"112 - m^^ - 8/ + 3/i^ - /i^ (8 - 8/i^ + /) C6 

32 - 152/1^ + 126/ - 34/i^ + 3/ - /(6 - 6/1^ + ^/) Cg 



+ 



/^^K,)' 



14 - IG/i^ + 3/i^(3 - ce)] log/i| 

^^^^^(l + 21og^)(l + 4m^C4) , 

48 rriN 654 
3/iM f 8-8/i^ + /i^ /iarccos(-|) / ^ ^x 2, 

(47rF^)H 4(4-/^2) (4-/i2)3/2 V '^ /- y r- &A- 

X (l + Cg + 2C7) , 



(3.26) 



where, again, 6^4 has the infinite part, which is given by 



A 



674 



3m Af ^ 



A-niNCi) , 



(3.27) 



(3.28) 



subtracted, while 654 in the isoscalar radius (in analogy to dj for the electric isoscalar radius) 
is finite. The isovector magnetic radius explodes like 1/M^ as the pion mass goes to zero, 
while, again, the isoscalar radius is finite in the chiral limit. These chiral limit results are 
to be expected because the Sachs form factors inherit the well known chiral singularities of 
the isovector Dirac and Pauli form factors, more precisely, the corresponding radii diverge 
in the chiral limit as 



„v\2 



[rlf 



lQg\ + 2 M 

log + 



(47rF^)2 



rriN 
+ ... , 



(3.29) 



where the ellipsis denotes subleading terms as the pion mass vanishes. This dominant pion 
loop effect in the nucleon form factors has already been established a long time ago [Q. It 
can also be expected on physical grounds: the pion cloud, which is Yukawa suppressed for 
finite pion mass, becomes long-ranged in the chiral limit leading to a divergent contribution. 
Such phenomena are also observed in two-photon observables like e.g. the electromagnetic 
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polarizabilities of the nucleoli. Also, as a consequence of analyticity, in the heavy baryon 
approach (i.e. in a strict chiral expansion) the isoscalar form factors are, to one-loop order, 
given simply by a polynomial subsuming tree and counterterm contributions. There are, 
however, isoscalar loop effects as higher order corrections, as is obvious from the above. 
These stem from the diagrams with the photon coupling to the nucleons (see graphs (5) 
and (10) in fig. ^, which only yield constant terms needed for renormalization in HBChPT, 
but also contain higher order t-dependent terms. As the spectral function corresponding to 
these diagrams only starts to contribute for t > 4r?2^, these t-dependent terms (in the real 
part) are highly suppressed. 



In section |3.2| , we already mentioned the marked difference in the threshold behavior of the 
isovector spectral functions in the relativistic and the heavy baryon approach. We only quote 
the results for the threshold behavior here and relegate the full results to appendix 0. Note 
that these have already been given in the literature, the relativistic results to third order 
in 101, the additional fourth order contribution as well as the heavy baryon expansions in 
[^ . We reemphasize that the imaginary parts in the infrared regularization approach do not 
differ from those in a fully relativistic calculation. For Im F^{t), the relativistic threshold 
expansion is given by 

and for ImF^(t), one finds 

ImF^it) = '^''!;^ it - AMlf/^ + Oiit - AMlfl^) . (3.31) 

In the heavy baryon scheme, however, the heavy mass expansion yields the expressions 



1 It - 4A/2)3/2 



which evidently violate the required threshold behavior ImFj^(t) ~ (t — 4M^)^/^ as pointed 
out in p. As remarked before, this is due to the fact that the two-pion threshold and 
the singularity on the second Riemann sheet, inherited from the vrvr -^ NN partial waves, 
coalesce in the heavy baryon limit. In contrast, the relativistic approach of course yields the 
correct analytical structures (for a more pedestrian discussion of this point, see e.g. [|T6[] ). 
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3.4 Results and discussion 

First, we must fix parameters. We use F.„ = 92.4 MeV, M.„ = 139.57 MeV, qa = 1.26, and 
m^v = {f^p + ^^n)/2 = 938.92 MeV. The LECs C2 and C4 are taken from the analysis of pion- 
nucleon scattering []I1|, |r^, C2 = 3.2 GeV^^ and C4 = 3.4 GeV~^. The LECs Cq and C7 can be 
pinned down from the well known magnetic moments of proton and neutron. The remaining 
two "electric" {d^, dj) and two "magnetic" LECs (654, 674) are determined from the electric 
and magnetic radii as given by the dispersion theoretical analysis of [Q, these values are 
0.847 fm, (rlY = -0.113 fm^, rL = 0.836 fm and rlr = 0.889 fm. We note that the 



r 



p 



E — >J-"'±i ^^'■h V E) — w.xxoiiii 1 I M — "J-o-JUiiii ciiiu / jvf 



long standing discrepancy of the proton charge radius determination from electron-proton 
scattering |^0| and Lamb shift measurements [^] has been resolved, converging to a value of 



rg = (0.88 ± 0.01) fm. We do not use this value here because so far no dispersion theoretical 
analysis exists using this novel information. 

As already stated, the Dirac and Pauli form factors are the natural quantities in any rela- 
tivistic approach like the one used here. However, for easier comparison to existing heavy 
baryon calculations and to the low energy data, which are often presented in terms of Sachs 
form factors, we will discuss these in detail. 

First we discuss the neutron charge form factor, as shown in fig. |^. The third and fourth order 
results compare well to the dispersion theoretical analysis and the trend of the novel data up 
to Q^ = 0.4 GeV^. The latter have been obtained using tensor polarized deuterium, polarized 



He and polarization transfer on deuterium |22|. Although it seems that the third order curve 



meets the novel data even better, we observe that the fourth order one is basically identical 
to the dispersion theoretical fit up to Q^ = 0.2 GeV^. Also given in that figure are the heavy 
baryon results to third 0] and fourth [^ order. These curves can easily be obtained from 



our analysis by performing the large mass expansion as explained before. Clearly, neither of 
these two curves is in acceptable agreement with the data, and, what is more, there is no 
sign of convergence so far. The resummation of the 1/m terms in the relativistic approach 
vastly improves the convergence of the chiral representation. This sensitivity of the neutron 



charge form factor to such recoil corrections was already anticipated in [^]. We remark 
however that the only fourth order contributions to the electric form factors in the heavy 
baryon approach are 1/m-corrections to third order diagrams, in other words, the full fourth 
order heavy baryon result for these can already be obtained from expanding the third order 
relativistic one. Therefore, the difference that does show up between the relativistic third 
and fourth order curves is really, in terms of strict chiral power counting, of higher order and 
hence (comparably) small. 

Consider the proton electric form factor next, shown in fig. ^. As in the neutron case, the 
third and fourth order curves are very close (where, however, the remarks made about the 
reasons for the smallness of this difference apply equally here), but here essentially show the 
linear behavior from the term proportional to the charge radius. Compared to the dispersion 
theoretical result (which describes the data very well in this range of momentum transfer), 
there is clearly not enough curvature. Polynomial terms of order t"^ (and higher) are not 
included up to fourth order, and the curvature coming solely from loop functions in the 



15 



one-loop approximation is obviously not sufficient. Fig. ^ also displays the results of the 
third p, |] and fourth [|23[ order heavy baryon calculation. While the general trend of the 



heavy baryon results is similar to what is obtained in the relativistic approach (far too small 
curvature to meet the data), we note that, while there is at least a slight improvement from 
the third to the fourth order relativistic calculation, the description visibly worsens for the 
fourth order heavy baryon result, as compared to the third order one. 

We now turn to the magnetic form factors of proton and neutron as shown in figs. ^, |[ The 
qualitative behavior of the different curves is very similar for both of them: To third order, 
the momentum dependence is given parameter-free because the LECs at this order only 
affect the normalization. We see that the 1/m corrections present in our approach visibly 
worsen the prediction obtained previously in the heavy baryon limit. That result is based on 
the leading chiral singularities given in eq. ( |3.29| ) because to this order, the corresponding 



isoscalar form factor is simply constant. We conclude that the recoil corrections are rather 
large and that the leading chiral limit behavior is not a good approximation for the Goldstone 
boson contribution to the magnetic radii in the case of finite pion masses. At fourth order, 
due to the presence of the counterterms from C^-lf, the radius can be fixed at its empirical 
value. Consequently, the full relativistic result and its heavy baryon limit are not very 
different any more. The absence of curvature terms of order t^ (and higher) is clearly visible 
in figs. §, §. 

The physics underlying this deficiency will be addressed in section ^. Only a uniformly 
reliable description for all four form factors is acceptable, and hence it is absolutely necessary 
to understand what is missing in the three dipole-shaped form factors: with large corrections 
for these, one might expect large corrections also for the neutron charge form factor, such 
that its very good description presented here might turn out to be only accidental. In the 
following, we shall show that this is indeed not the case. 

4 Inclusion of vector mesons 

It is well established that the vector mesons contribute significantly to the nucleon form 
factors. For example, extended unitarity allows one to reconstruct the isovector spectral 
functions |^ below t ~ 1 GeV^ from the pion form factor and the analytically continued 



TTTT -^ NN isovector partial wave amplitudes. Besides the important contribution of the 
two-pion continuum, the p meson clearly shows up. Similarly, in the isoscalar channel, the 
uj and the dominate the spectral function at low positive t. All these effects are clearly 
visible in dispersion theoretical analyses of the nucleon form factors, see [^ |, ^. In the 
chiral expansion performed in the preceding sections, such effects are included in the low 
energy constants. This follows directly from the low momentum expansion of any vector 
meson propagator (for simplicity, we do not show the explicit Lorentz structures). 
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In a fourth order chiral analysis as presented here, one is sensitive up to the terms hnear 
in t, which contribute to the various electromagnetic radii. Stated differently, any vector 
meson contribution is hidden in the fit values of the various LECs. As we discussed before, 
the curvature effects on the electric and magnetic form factor of the proton as well as the 
magnetic neutron form factor are much too small. This can be cured in two ways. One could 
either attempt a higher order analysis or include vector mesons dynamically. While the first 
approach would be more systematic, we choose here the second, for various reasons. With 
explicit vector mesons (VM), we not only account for all terms in the expansion eq. ( [4.1[ ) 
but also do not introduce any new unknown parameter. This can be understood as follows: 
adding the vector mesons in a chirally symmetric manner and retaining the corresponding 
dimension two, three and four counterterms, any LEG «« takes the form 

a, ^ Aj + vector meson — contribution , (4-2) 

where the remainder dj parameterizes the physics not related to the explicitly included vector 
meson contribution. For such a decomposition to make sense, the scale dependent LECs 
should be taken at A = My- Since in the infrared regularization procedure, the nucleon 



mass is taken to be the intrinsic scale for all loop integrals as described in section ^]2 
here we neglect scale mismatch due to My 7^ rriN (which is probably justified because of 
log{mM/Mp) ^ 0.2). The vector meson propagator generates a whole string of higher order 
terms, hopefully resumming the most important contributions not included at fourth order. 
The parameters appearing in the vector meson contributions will be taken from existing 
dispersive analyses of the nucleon electromagnetic form factors. We therefore only need 
to refit the low energy constants Aj. This also allows to study the concept of resonance 
saturation. It was already shown in |^6[ that the numerical values of the dimension two 
LECs can be understood from s-channel baryon and t-channel meson resonance excitations, 
in particular from the A(1232) and the p. In the case at hand, we can investigate this 
concept for the LECs Cg, C7, d^, d^, 654, and 674. It was already noted in |2^ that cg and C7 
are largely saturated by vector meson contributions. 

4.1 Chiral Lagrangians including vector mesons 



We employ the tensor field representation of spin-1 fields as advocated in |^, i.e. the vector 
mesons are written in terms of antisymmetric tensor fields 

W^, = -W,, (4.3) 

with the three degrees of freedom Wij {i,j = 1, 2, 3) frozen out. This representation is most 
natural for constructing chirally invariant couplings of vector mesons to pions and photons 
because no particular dynamical character of the vector mesons is assumed (in contrast e.g. 



to the massive Yang-Mills or hidden symmetry approaches as reviewed in p8| , |29|.) In this 
section, we temporarily switch to SU(3) chiral Lagrangians, because the (p meson contributes 
to the isoscalar form factors. The free Lagrangian then takes the form 

1 M^ 

Cy = — d^wn^OpW""^^ + — ^iy;^iy^^''^ , (4.4) 
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where 



W, 



flV 



V K*- 






K*+ \ 

-0 / 



(4.5) 



^u 



Here, we have written the vector meson matrix in terms of physical particles, assuming ideal 



'-UJ mixing. From eq. ( [4.4|) , one easily derives the vector meson propagator as given in |^ 



G 



liu,pa 



'x,y) 



{0\T{W,,{x),W,^{y)}\0) 



d^k 



Jk-{x-y) 



M^J (27r)4 M^ - P - ie 



(4.6) 



According to [^ , the lowest order interaction with Goldstone boson fields as well as external 
vector and axial-vector sources can be written as 



1 



^w = ^[FviW^'^p;^) + zGviW^'-l 



(4.7) 



where the couplings Gy and Fy can e.g. be determined from the decay widths p — > tttt and 
p — > e+e~. Following [^, the lowest order couplings of massive spin-1 fields to baryons can 
be written in terms of a chirally invariant Lagrangian as 

+Sn/F {BrUD", W,,],B)±) + Ss {B^B) ([D^ W^,]) 

+Ud/f {Ba"'iW,,,[Dx,[D'',B]])^) + Us {Ba'^lDxAD^, B]]){W,,) , (4.8) 

where, as usually in SU(3), the index D refers to the anticommutator {A,B)+ = {A,B}, 
while the index F accompanies the commutator [A, B)_ = [A, B]. In addition, we introduced 
singlet couplings with indices S. These coupling constants are related to the ones used in 
the Lagrangian of the standard vector representation of the p, 



^pN= -gpNN'^{Ypp- 



2m 



V^^/.. 



N 



duPp-r}^ 



by 



QpNN 



QpNN l^p 



^{mN{UD + Up)-2iSn + Sp: 



M. 



In analogy, one finds for the cj and (p couplings 

My/ 



9uiNN 



-^{niNiUD + Uf + 2Us) - 2{Sd + Sf + 23. 
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(4.9) 

(4.10) 
(4.11) 

(4.12) 



g.NNK^ = -^^^^{Rd + Rf + 2Rs) , (4.13) 

g^NN = -My{m^{UD -Uf + Us) - 2{Sd -Sf + Ss)) , (4.14) 

g4>NN K.4, = -n—i^D - Rf + Rs) ■ (4.15) 

Mv 

In PD|, one additional term (proportional to new couplings T^/p) of higher order in deriva- 
tives was introduced, 

-^W = Td/f {Br{\Dx, W^X [D\ \D\ S]])±) + Ts {Br[D\ [D^ B]]) {[D^, W^,]) , 

(4.16) 
which was subsequently shown to contribute to the electric form factors as a term 0{q'^), 
therefore beyond the order to which we are working here; the authors of ||30| however left 
out a term 

'^Ibw = Vd/f {Ba^\[Dx. [D^\ W^X,B)±) + Vs {Ba^^B) ([Da, [Z^^ W,,]]) (4-17) 

which enters the magnetic radius, i.e. contributes a term 0{q^) to the magnetic form factors 
and hence should be considered in a 0{q'^) calculation. However, as nothing is known from 
elsewhere about such couplings, we shall not introduce these here.p^ Note that, via the 



g^-expansion of the resonance pole, the contributions to the magnetic moments stemming 
from the Rd/f/s terms do of course induce additional g^-dependence in the magnetic form 
factors. 

There is no generalization of chiral perturbation theory which fully includes the effects of 
vector mesons as intermediate states to arbitrary loop orders. As the masses of the vector 
mesons do not vanish in the chiral limit, they introduce a new mass scale which, when 
appearing inside loop integrals, potentially spoils chiral power counting in a similar manner 
as the nucleon mass in a "naive" relativistic baryon ChPT approach. In analogy to the heavy 
fermion theories, "heavy meson effective theory" has been used to investigate vector meson 
properties like masses and decay constants when coupling these to Goldstone bosons in a 
chirally invariant fashion (see [Q or for some special processes [^). This approach does 
not work with the heavy particle number not conserved as in processes with these resonances 
being virtual intermediate states. However, these difficulties do not arise as long as there is 
no loop integration over intermediate vector meson momenta. Indeed, to the order we are 
working here, we can even set up a "power counting scheme" for diagrams including vector 
mesons which allows to calculate corrections to the simple tree diagrams (where the photon 
couples to the nucleon via an intermediate vector meson). We count 

1. the couplings of vector mesons to the photon and of the p to two pions as O(g^) (see 
eq. ( [4.7|) and the usual power counting for F^ and u^); 

2. the tensor coupling of vector mesons to nucleons {Rd/f/s in our notation) as 0{q^), 
the vector coupling {Sd/f/s, Ud/f/s) as 0{q^); 



'^^The terms in eq. ( |4.16| ) were allowed for in the analysis |Q of TriV-scattering and found to be very small. 
This lends credit to the omission of higher order couplings here. 
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3. the vector meson propagator as C(g°) (see eq. ([4.6D). 

Up to 0(5'^), we therefore have to include the diagrams shown in fig. |^. The numbering 
indicates the correspondence to diagrams with contact terms as shown in fig. ^j: the tree 
level coupling of the vector mesons to the nucleon via the tensor coupling, diagram (2*), 
enters the magnetic moments as do the LECs cq, cj in diagram (2), diagrams (10) and (10*) 
as well as (11) and (11*) yield vertex corrections to these due to pion loops, the vector 
coupling in diagram (3*) compares to the LECs (ig, d-j as shown in diagram (3), and finally 
the p-exchange in diagram (12*) contributes to the vrvrA^A^-coupling ~ C4 in diagram (12). 

Explicitly, the analytic results for the form factors including vector mesons can be gained 
from those with contact terms only by the following replacements: 

C6 -> h + 9pNN K.p ^2 _^ ' (^•^^) 

- _ 9pNN i^p FpMp Qujnn i^uj F^M^^ g(t,NN 1^4, FfffM^ 
ci -^ C7 ^ mF^ 2 'W^t^ 2 mJ^' ^ ^ 

p U) (p 

„■ -jr 9pNN Fp 1 

''^ "^^ 2 Mp M^^-t' ^^-^^^ 

Especially, these hold for the loop diagrams with pion loops as vertex corrections to photon 
couplings via Cq, Cj, see diagrams (10), (11) in fig. || and (10*), (11*) in fig. ^. At leading 
order, resonance saturation for cq, c-j has already been investigated in |^. The vector 
meson contributions to the magnetic moments and the electric radii can be found trivially 
from eqs. ( |3.19| )-( p.24| ) by using eqs. ( [4.18| )- p.21| ) at t = 0. The vector meson contributions 



to the magnetic radii, finally, yield a more complicated replacement law for 654, 674 due to 
the aforementioned loop corrections, which, however, at leading order reads 

If F^ Fs^ 

654 -^ ^54 + ^ |fi'c^AfiVft:a;T^ + 5'0JViV/^<ATT3| + C'(/i^) , (4.22) 

674 ^ e'^^ + ^^9pNNKpJ^ + 0{^^^) . (4.23) 



Finally, resonance saturation for the LEG C4 has also been analyzed in |^, where it was 
found that C4 is completely saturated by p. A, and (a very small) Roper contribution. Here, 
we only want to replace the p contribution by its dynamical analogue, therefore setting 



2 



,,^e. + M£^^A-8. + ^.#^. (4.24) 



2mN M^-t AruN M^ - t 



where, in the second step, a universal p-hadron coupling g = gpNN = gpmr = GpMp/F^ and 
the KSFR relation M^ = 2g'^F^ [^ were assumed. 
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4.2 Results and discussion 

Before presenting results for the four form factors, including the effects of dynamical vector 
mesons, we have to discuss the values for the miscellaneous coupling constants introduced 
thereby. First of all, for the vector meson masses we use Mp = 770 MeV, M^ = 780 MeV, 
M^ = 1020 MeV. The couplings to the photon have been fixed from the partial widths of 
V -^ e+e- (in analogy to §) to be Fp = 152.5 MeV, F^ = 45.7 MeV, F^ = 79.0 MeV, 
the ratios being in satisfactory agreement with what one would expect from exact SU(3)- 
symmetry (plus ideal (p-LU mixing), Fy = Fp = 3F^ = S/v^F^. 

The couplings of vector mesons to nucleons are taken from the most recent dispersive analysis 
[|]. The respective values (adjusted to our conventions) are 



(4.25) 



The resulting electric and magnetic form factors of the proton and the neutron are shown in 
figs. ^HTTl Consider the electric form factors first. The vector meson contribution already 
supplies sufficient curvature for an adequate description of the proton charge form factor 
at third order, cf. fig. ||. This is achieved without new adjustable parameters, but simply 
due to the higher order terms induced by the inclusion of the vector mesons. It is worth 
noting that at fourth order, the chiral plus vector meson representation is in almost perfect 
agreement with the result of the dispersive analysis up to momenta of Q^ = 0.4 GeV^, i.e. 
for much higher momenta than considered so far in chiral perturbation theory approaches. 
However the already well-described neutron electric form factor, see fig. ^, is, at fourth 
order, only very mildly affected by the vector meson contribution, but is of course also most 
sensitive to the precise choice of the vector meson couplings. We only want to state that for 
a reasonable choice of these parameters, the good result of the chiral one-loop representation 
is not spoiled.p^ 



dpNN = 


4.0, 


Kp = 


6.1 , 


QljNN = 


41.8 , 


I^UJ ^ 


-0.16 , 


g<l>NN = 


-18.3 , 


f^<l, = 


-0.22 . 



Turning now to the magnetic form factors, see figs. |T0|, |TT|, we find that both the third and 
the fourth order curves yield reasonable descriptions of the data, in both cases the fourth 
order results being slightly better than the respective third order ones. The latter show 
the right slope at the origin (which is fitted in the fourth order case), and all curves have 
approximately the right curvature, which is a bit stronger at fourth order. 

We have, in addition, plotted the third and fourth order curves for G^, G^^, and G^ 
(including vector mesons), normalized and divided by the dipole form factor 

Gd{Q') = - ^ -2 , (4.26) 

1 + gV0.71 GeV^ 



*^ We also note that the dispersive analysis of ^ does not include most of the new data Q . Refitting the 
vector meson parameters accordingly to the new data base will also lead to an improvement of the fourth 
order curve in our analysis. 
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see figs. p!^-p^. Here, again, we see that the fourth order prediction for the proton charge 
form factor agrees extremely well with the data up to Q^ = 0.4 GeV^. For the neutron 
magnetic form factor, the fourth order curve, though rising above dipole and dispersion 
theoretical fit, still is within the error range given by the data, while the fourth order result 
for the proton magnetic form factor deviates from the data above Q^ = 0.3 GeV^. 

It is worth pointing out the following essential conclusions drawn from this discussion: 



P \2 _ 

m) - 


- vmJtt + (^m)vm 


n \2 
M) 


= vm)tt + (^m)vm 



1. In section |3.4| we saw that the "prediction" of the magnetic radii in a third order cal- 
culation worsens considerably when comparing the relativistic and the heavy baryon 
approach. This observation indicated that the large contribution of the pion cloud 
to the magnetic radii was an artifact of the 1/m-expansion to 0{q^). Inspection of 
the third order curves including explicit vector meson effects shows that this reduc- 
tion of the pion cloud effect is indeed consistent with the fairly conventional vector 
meson parameters, yielding in total excellent results for the magnetic radii which are 
still "predicted" and not fitted. We here show these "predictions" with the separate 
contributions of pion-cloud and vector mesons, 

(0.18 + 0.49) fm^ = (0.82 fm)2 , (4.27) 

(0.26 + 0.44) fm^ = (0.84 fm)^ , (4.28) 

in good agreement with the values from dispersion analyses, r^ = 0.836 fm and 
r^ = 0.889 fm. Note however that, although we did not fit any new parameters 
to the magnetic radii, additional experimental information enters via the vector meson 
couplings which in turn have been fitted to the form factors in the dispersive approach. 

2. This analysis indicates that a "strict" ChPT calculation to higher-than-fourth order is 
probably only moderately sensible. Apart from various two-loop contributions which 
amount only to vertex corrections of diagrams already present in the one-loop case, 
the major "new" contribution is the three-pion continuum entering the isoscalar form 
factors, which was already shown to yield a negligibly small correction to the uj/(f)- 
peak in the isoscalar spectral function [Q. Therefore, an analysis to 0{q^) or 0{q^) 
would basically result in fitting new counterterm parameters in order to reproduce 
the resonance contributions considered explicitly here. Apart from that, it is so far 
not known how to generalize the infrared regularization scheme to higher loop orders, 
and therefore for such a higher order calculation one would have to retrieve to the 
non-relativistic formalism. 

4.3 Remarks on resonance saturation 

As we have not assumed resonance saturation of the various low energy constants, but taken 
the resonance parameters from elsewhere and refitted the contact terms, we may now check 
how well resonance saturation by the lowest lying vector mesons actually works. Table |l| 
shows the values for the low energy constants, fitted at third and fourth order, with and 
without vector meson contributions. The vector meson contribution to the LEG C4 as given 
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<l' 


q' 


q^+VM 


g^+VM 


c& 


5.03 


4.77 


0.20 


-0.06 


Cl 


-2.72 


-2.55 


-0.27 


-0.09 


dl 


0.67 


0.74 


1.34 


1.40 


dj 


-0.70 


-0.69 


-0.04 


-0.03 


654 


— 


0.26 


— 


-0.04 


674 


— 


1.65 


— 


-1.15 



Table 1: Values for the various low energy constants, with and without vector meson con- 
tributions. Note that the third and the fourth columns, strictly speaking, refer to the LECs 
dj instead of aj. Cg, C7 are dimensionless, the di are given in GeV~^, the Cj in GeV~^. 



in eq. (f4.24|) numerically amounts to 1.6 GeV~^ (see also ||26|) which reduces the value of C4 
by about 50%. This remainder which is due to (mainly) A and a small Roper contribution 
is still represented by a contact term. 

We conclude that the low energy constants entering the magnetic moments (cg, C7) as well as 
those fitted to the isoscalar radii, both electric (^7) and magnetic (654), are nearly perfectly 
saturated by the lowest lying vector meson nonet, while the contact terms entering the 
isovector radii {d^, 674) are not at all. This can be interpreted to the effect that u and 
are already sufficient to describe the isoscalar channel of the vector form factors to good 
accuracy, while higher resonances are mandatory for an adequate description of the isovector 
form factors. This is in agreement with what is found in dispersive analyses. However this 
does not necessarily invalidate our approach, as these higher resonances have much larger 
masses (the p(1450) being the lowest lying of them) and therefore will hardly contribute to 
the curvature of the form factors (which is the change we are looking for, compared to linear 
contact terms). 



5 Summary 

We have studied the electromagnetic form factors of the nucleon in a manifestly Lorentz 
invariant form of baryon chiral perturbation theory to one-loop (fourth) order. As discussed, 
in this scheme based on the so-called infrared regularization of loop graphs, one is able to 
set up a systematic power counting scheme in harmony with the strictures from analyticity. 
The pertinent results of our investigation can be summarized as follows: 

(1) To fourth order, the neutron and proton electric form factors each contain one low- 
energy constant which can be fixed from the empirical information on the corresponding 
radii. This gives a good description of the neutron charge form factor up to four- 
momentum transfer squared of Q^ = 0.4 GeV^ and, furthermore, exhibits convergence 
in that the corrections when going from third to fourth order are small. This is in 
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contrast to the heavy baryon expansion and can be traced back to the proper resum- 
mation of the recoil terms in the relativistic expansion. For the electric form factor 
of the proton, the one-loop representation gives too little curvature and thus deviates 
from the data already at Q^ ~ 0.2 GeV^, similar to the heavy baryon description. 
However, no large fourth order corrections are found below Q^ = 0.4 GeV^. 

(2) To third order, the momentum dependence of the magnetic proton and neutron form 
factor is given parameter-free. The 1/m corrections present in our approach worsen 
the prediction for the magnetic radii based on the leading chiral singularities, like 
e.g. in the heavy baryon approach. The leading chiral limit behavior is not a good 
approximation for the Goldstone boson contribution to the magnetic radii. At fourth 
order, the magnetic radii can be fixed. Again, there is not enough curvature in the 
one-loop representation and one observes large corrections when going from third to 
fourth order already at Q^ — 0.1 GeV^. 

(3) We have demonstrated explicitly that the spectral functions of the isovector form 
factors have the correct threshold behavior. The strong momentum-dependence of 
these spectral functions close to threshold is due to the branch point singularity on the 
second Riemann sheet inherited from the tttt -^ NN P-wave partial wave amplitudes. 

(4) We have included the low-lying vector mesons p, uj, (p in a chirally symmetric manner 
based on an antisymmetric tensor field representation. This does not introduce any new 
parameters since these (masses and coupling constants) are taken from the PDG tables 
and from a dispersion theoretical analysis. Refitting the previously defined low-energy 
constants by subtracting the vector meson contribution, we find a good description 
of all four form factors already at third order, with small fourth order contributions, 
which further improve the theoretical description. In particular, we demonstrate that 
the vector meson contributions cancel to a large extent in the neutron charge form 
factor, thus solidifying the result obtained in the chiral expansion. 

(5) The inclusion of vector mesons allows to investigate the resonance saturation hypothesis 
for these couplings. We find that the couplings related to the magnetic moments and 
the isoscalar radii are almost completely saturated by the low-lying vector mesons. 
This is, however, not the case for the LECs entering the isovector radii. This can be 
traced back to the fact that while the u and the already give a good description of 
the isoscalar form factors, for the isovector ones one has to include higher mass states 
than the p, in agreement with findings from dispersion theory. 
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A Loop integrals 

In this appendix, we define the loop integrals needed in this paper and evaluate them in the 
infrared regularization scheme. Several of these results have already been given in 0. 

A.l Definition of the loop integrals 

We use the following notation: 

M^ t . t 

niN rn% M^ 

In addition, everywhere except in the loop integrals with just one meson and one nucleon 
propagator, we only need the case where the nucleon momenta are on-shell, i.e. p^ = p'^ = 



Define the following loop integrals: 

d'^k 1 



i J I (27r)'^ Ml - F 
d'^k 1 



A. , (A.l) 



i J I {2tiY [Ml - A;2] [Ml - {k + q)^] 

1 r d'^k k^ _ 1 

i J I (27r)^ [Ml - P] [Ml -{k + g)2] ~ ~2 



At) , (A.2) 

q^J{t) , (A.3) 



I J I {2^Y [Ml - k^][Ml -{k + g)2] 

1 r d'^k 1 

7 J I {27iY [Ml - A;2] [m% - {p - k^] 

1 r d'^k kf, 

1 J I {2ttY [Ml - fc2] [m% -{p- kf] 



I r d^k 




1 






I J I (27r)'^ [Ml - 
I r d'^k 


- fc2] [m\i - 


- {p - kY][mlf - 


-ip'- 


-kY] 


I J I {2t^Y [Ml - 
I /• d'^k 


- A;2] [77^2^ - 


- {p - kY][m% - 


-ip'- 


-kY] 



i J I {27iY [Ml - fc2] [m% -ip- kY] [m% - {p' - kf 



%(iu - Qi^ut) A^^ (t) 

+ q,q.A'\t) , (A.4) 

= /(/) , (A.5) 

= P,I^'\p') , (A.6) 

= Ut) , (A.7) 

= Q,I^2\t) , (A.8) 

+ Q,Qulf{t) 

+ q.qulfit) , (A.9) 
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d'^k 1 



i J I (27r)^ [M2 - A;2] [M^ - {k + q^] [m% - {p - k)^] 
1 r d'^k k„ 



i J I (27r)^ [M2 - A;2] [M^ - (A; + g)^] [m^ - (p - A;)2] 



/2i(t) , (A.IO) 

^g^J2i(t) , (A.ll) 



1 f d k k^ky ^ QQ 



i J I {2nY [Ml - A;2] [M^ - {k + qY] [m% - {p - k) 
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9,.iTAt) 



QQ/ 



+ q.q.iim (A.12) 

where Jj symbolizes loop integration according to the infrared regularization scheme. 

A. 2 Reduction of the tensorial loop integrals 

The reduction of the tensorial loop integrals to the corresponding scalar ones can be per- 
formed in the standard way and leads to the following results: 

j'^'Ht) = ^^^l^{(t-4M,^)J(t) + 2A.} , (A.13) 

J(2)(t) = \jit)-Yt^. , (A.14) 

I^'\p') = ^{(p'-< + MMp')+A.} , (A.15) 

I^2\t) = ^-^^{limD + M^lAit)} , (A.16) 

I^A\t) = j^{Ut) - I^A\t)}M'^ , (A.17) 

I^I\t) = (^ _ ^^(^^^ _ ^^ {{{d - l)li\t) - lA{t))Ml + ^/(^)(0} , (A-18) 

ifit) = JJ^t{{li\t)-lAit))M'^-^I^'\ml)] , (A.19) 

I?i{t) = ^^^^^^^-^{(2M2-t)/2i(t)-2/(m^) + 2J(t)} , (A.20) 

/2TW = —^{2I{ml)-{4M'^-t)h,{t) + 2{2M'^-t)lg{t)} , (A.21) 



4(2 - d) 



'^^ = 4(.-2)(U-i) {"""'^'-^"'-^>^"'' 



rrir 



-ml - t)hi{t) + 2{d - l){2Ml - t)lg{t) \ , (A.22) 
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Iim = ^^^l^{-2(rf-3)/(m^)+2(d-2)/«(m^) 

-(4M2-(rf-l)t)/2i(t) + 2(2M2-t)/2«i(t)| . (A.23) 

A. 3 Scalar loop integrals 

The scalar loop integrals are found to be 

A. = 2M2|L + ^log/i} , (A.24) 



J{t) = -2{^ + I^log4-I^(l + Mt)) , (A.25) 

Iim%) = V(^ + Y^log/i) + Y|^{|-y4^arccos(-0} , (A.26) 



where the following loop functions have been reduced to integrals over one Feynman param- 
eter: 



k{t) = /^^xlog(l-x(l-x)g) = y^^log( ^^+^ )-2 , (A.29) 

/W = / ^a;- ^= . ^. ^logp^ ^ ^ , A.30 

Jo 1 — a;(l — x) r V" ^(4 — t) Vv4 — r — v— t"/ 



-1 arccos ^ 



^(t) = / dx- , ; , \ , (A.31) 

/•I logfl — x(l — x)9) 
hit) = / dx^ -^ -^^ , (A.32) 

Jo 1 — X(l — X)T 

, arccosf- , ^(^-(^-)^) ) 

/..(t) = / dx- WoEH^^^SL . (A.33) 

•^0 (l-x(l-x)r)^l-^-x(l-x)^ 

B Form factor contributions from separate diagrams 

In this section, we give the contributions to the form factors Fi{t), F2(t), coming from the 
various diagrams shown in fig. || . 
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B.l Contributions to Fi 

Ei+3 = l-{r'^d(i + 2d-r)t , (B.l) 

+8m%I^]\t) + 32m%I^^\t)} , (B.2) 

E, = -^^r'[tj(^)(t)+4m%I^',{t) + 16m%lg'^{t)} , (B.3) 

n 

E7 = ^r'{K-2mll^'\m%)} , (B.4) 



Es = -^^^ > (B.5) 

E, = ^,tJ^'\t) , (B.6) 

n 

E,o = ^((3-r3)c6 + 6c,)t/Jf)(t) , (B.7) 



B.2 Contributions to F, 



2 



M2+3+4 = - (1 + r^) C6 + C7 + (r^ 4 + 2 ^7) t + 2m^ (2 654 + r^ 674) t 

-8m^M2(2eio5 + r3eio6) , (B.9) 



M5 = -||(3-r^)4m^^/J)(t) , (B.IO) 

Me = ||r3l6m^/g«(t) , (B.ll) 

IT 

Mio = -^((3-r3)c6 + 6c7){A.-4my«(m^)+4m^M2j^(t) 

-16m^/f(t) + 8m^t(/(f)(t)-/(f)(t))} , (B.12) 

M,2 = l^r»C4iJ<'>(«) . (B.14) 



7r 
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B.3 Z-factor 

We here spell out the Z-factor needed for wave function renormalization up to fourth order. 
For the evaluation of the form factor Fi, the nucleon charge has to be multiplied by Z^v, 
while for F2, the anomalous magnetic moment (a second order contribution) only has to be 
renormalized by the Z-factor up to third order, i.e. the term ~ C2 can be dropped. 

, ^3M,^ / Ml^ 

C Imaginary parts and spectral functions 

Out of the diagrams depicted in fig. |, only graphs (6), (9), and (12) contribute to the 
spectral functions of the isovector form factors in the low energy region, i.e. starting at the 
threshold t = 4M^. The separate contributions can be calculated either from the imaginary 
parts of the two basic loop functions involved (see also |^), 



1 / 4M2 



ImJ2i(t) = = arctan —-: — , (C.2) 

Svr ,,/t(4m^ - t) t-2M^ 



or directly by using Cutkosky rules. This leads to the following expressions: 
T^gl r 4M^ 



Imi^fi 



X 



1927rF2 {Am% - tf \ t 

I lQm% {ht - SMl) + Ami t (st - 14M^) -t^{t- AMl) (C.3) 



mU?>e - l2Mlt + 8M4 + Mtt J {Ami - t){t - AM^) ] 



47rF2 {Ami - i)M ^ I ^ 



(2ml + t){t-AM^)+6M^ .^{Aml - t){t - AM^^)\ 

+2 ^ ^ ^ ^ arctan -^ —: } , C.4) 

/{Aml-t){t-AM^) t-2Ml J 



r^ it - 4M2?/2 
1927rF2 ti/2 ' ^ > 

For these results, see also §, |^ . 
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(b) 



Figure 1: Lowest order loop diagrams contributing to the mass renormalization of (a) the 
pion at 0{(f^)^ and (b) the nucleon at Oic^y Solid and dashed lines refer to nucleons and 
pions, respectively. 
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Figure 2: Feynman diagrams contributing to the electromagnetic form factors up to fourth 
order. Solid, dashed, and wiggly lines refer to nucleons, pions, and the vector source, re- 
spectively. Vertices denoted by a heavy dot / a square / a diamond refer to insertions from 
the second / third / fourth order chiral Lagrangian, respectively. Diagrams contributing via 
wave function renormalization only are not shown. 
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Figure 3: The neutron electric form factor in relativistic baryon chiral perturbation theory 
(sohd hues) to third (blue curve) and fourth (red curve) order. For comparison, the results 
of the heavy baryon approach are also shown (blue/red dashed line: third/fourth order). 
All LECs are determined by a fit to the neutron charge radius measured in neutron-atom 
scattering. Also given is the result of the dispersion theoretical analysis (black dot-dashed 
curve). The data are from [p2 |. 
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Figure 4: The proton electric form factor in relativistic baryon chiral perturbation theory 
(sohd hues) to third (blue curve) and fourth (red curve) order. For comparison, the results 
of the heavy baryon approach are also shown (blue/red dashed line: third/fourth order). 
All LECs are determined by a fit to the proton charge radius as given by the dispersion 
theoretical result (black dot-dashed curve). 
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Figure 5: The proton magnetic form factor in relativistic baryon chiral perturbation theory 
(sohd hues) to third (blue curve) and fourth (red curve) order. For comparison, the results of 
the heavy baryon approach are also shown (blue/red dashed line: third/fourth order). The 
third (fourth) order LECs are determined by a fit to the proton magnetic moment (radius). 
Also given is the dispersion theoretical result (black dot-dashed curve). 
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Figure 6: The neutron magnetic form factor in relativistic baryon chiral perturbation theory 
(sohd hues) to third (blue curve) and fourth (red curve) order. For comparison, the results of 
the heavy baryon approach are also shown (blue/red dashed line: third/fourth order). The 
third (fourth) order LECs are determined by a fit to the neutron magnetic moment (radius). 
Also given is the dispersion theoretical result (black dot-dashed curve). 
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Figure 7: Feynman diagrams including explicit vector meson contributions to the electro- 
magnetic form factors up to fourth order. Solid, dashed, double-dashed, and wiggly lines 
refer to nucleons, pions, vector mesons, and the vector source, respectively. The vertex de- 
noted by an open dot refers to the vector coupling of the vector mesons to the nucleon which 
is of subleading chiral order as compared to the tensor coupling. 
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Figure 8: The proton electric form factor in relativistic baryon chiral perturbation theory 
including vector mesons (solid lines) to third (blue curve) and fourth (red curve) order. 
For comparison, the results without vector mesons are also shown (blue/red dashed line: 
third/fourth order). All LECs are determined by a fit to the proton charge radius as given 
by the dispersion theoretical result (black dot-dashed curve). 
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Figure 9: The neutron electric form factor in relativistic baryon chiral perturbation theory 
including vector mesons (solid lines) to third (blue curve) and fourth (red curve) order. 
For comparison, the results without vector mesons are also shown (blue/red dashed line: 
third/fourth order). All LECs are determined by a fit to the neutron charge radius. Also 
given is the result of the dispersion theoretical analysis (black dot-dashed curve). The data 



are from [22 
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Figure 10: The proton magnetic form factor in relativistic baryon chiral perturbation theory 
including vector mesons (sohd hues) to third (blue curve) and fourth (red curve) order. 
For comparison, the results without vector mesons are also shown (blue/red dashed line: 
third/fourth order). All LECs at third (fourth) order are determined by a fit to the proton 
magnetic moment (radius). Also given is the dispersion theoretical result (black dot-dashed 



curve . 
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Figure 11: The neutron magnetic form factor in relativistic baryon chiral perturbation theory 
including vector mesons (sohd hues) to third (blue curve) and fourth (red curve) order. 
For comparison, the results without vector mesons are also shown (blue/red dashed line: 
third/fourth order). All LECs at third (fourth) order are determined by a fit to the neutron 
magnetic moment (radius). Also given is the dispersion theoretical result (black dot-dashed 



curve . 
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Figure 12: The proton electric form factor in relativistic baryon chiral perturbation theory 
including vector mesons to third (blue curve) and fourth (red curve) order, divided by the 
dipole form factor. For comparison, we show the dispersion theoretical result (black dot- 
dashed curve) and the world data available in this energy range. 
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Figure 13: The proton magnetic form factor in relativistic baryon chiral perturbation theory 
including vector mesons to third (blue curve) and fourth (red curve) order, divided by the 
dipole form factor. For comparison, we show the dispersion theoretical result (black dot- 
dashed curve) and the world data available in this energy range. 
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Figure 14: The neutron magnetic form factor in relativistic baryon chiral perturbation theory 
including vector mesons to third (blue curve) and fourth (red curve) order, divided by the 
dipole form factor. For comparison, we show the dispersion theoretical result (black dot- 
dashed curve) and the world data available in this energy range, where the data points 
denoted by squares (instead of circles) refer to the more recent measurements [0. The 
older data can be traced back from [Bl. 
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